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ABSTRACT 


This  paper  derives  an  algorithm  for  use  by  an  airborne 
mid-air  collision  avoidance  system  to  determine  when  an  alarm 
should  be  given  in  case  a mid-air  collision  is  imminent.  The 
algorithm  is  based  on  an  extension  of  the  standard  modified  tau 
alarm  criterion  used  in  most  collision  avoidance  system  threat 
logics. 

The  standard  criterion  uses  only  altitude  and  range  data 
and,  as  a result,  will  generate  high  alarm  rates  in  heavy  air 
traffic.  The  criterion  presented  here  makes  use  of  bearing 
data  as  well  as  altitude  and  range  data  and  should,  therefore, 
provide  lower  alarm  rates. 
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I.  INTRODUCTION 

The  objective  of  this  paper  is  to  construct  a new  algo- 
rithm for  deciding  when  the  danger  of  collision  between  two 
aircraft  is  great  enough  to  warrant  an  avoidance  maneuver.  In 
making  such  a decision  a collision  avoidance  system  (CAS)  in 
one  of  the  aircraft  would  measure  the  other's  relative  range, 
altitude,  and  bearing*  to  which  data  it  would  then  apply  the 
algorithm.  If  the  result  were  consistent  with  the  possibility 
of  an  imminent  collision  the  CAS  would  generate  an  alarm  to 
alert  the  pilot. 

The  scope  of  the  work  reported  here  has  been  limited  delib- 
erately to  the  mathematical  development  and  justification  of 
the  collision  threat  algorithm.  No  attempt  has  been  made  to 
assess  its  potential  usefulness  in  practice,  e.g.,  for  reducing 
the  number  of  false  alarms  relative  to  the  number  generated  by 
other  threat  criteria  which  have  been  implemented  or  suggested. 
Hopefully,  this  can  be  done  in  a future  study. 

ANTC-117  (Ref.  1),  the  document  containing  the  CAS  perform- 
ance standards  recommended  by  the  Air  Transport  Association, 
specifies  alarm  criteria  which  do  not  involve  bearing  data. 

These  criteria  are  of  two  types:  one  is  concerned  entirely 

with  altitude  and  altitude  rate  while  the  other  is  concerned 
entirely  with  range  and  range  rate.  The  second  is  based  on  the 
so  called  modified  tau  condition  which.  Ref.  2 has  shown, 

*The  relative  range,  altitude  and  bearing  rates  are  also  needed, 
but  these  quantities  can  be  calculated  from  measurements  of 
the  corresponding  static  quantities  at  two  different  times 
separated  by  an  accurately  determined  time  interval. 
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I provides  the  most  efficient*  possible  alarm  criterion  using 

range  and  range  rate  data  alone. 

However,  this  result  of  Ref.  2 is  valid  only  when  the 
fixed  parameters  of  the  modified  tau  condition  are  properly 

• selected.  Some  IDA  studies  (Ref.  3 and  Ref.  4)  have  applied 
the  analysis  of  Ref.  2 to  show  that  the  fixed  parameters  se- 
lected for  the  alarm  criteria  specified  in  ANTC-117  have  values 
which  will  not  guarantee  an  alarm  for  every  possible  threat. 

• In  addition,  the  studies  concluded  that  in  high  density  traf- 
fic such  as  projected  by  the  FAA  for  the  Los  Angeles  Basin  in 
1982  use  of  the  modified  tau  condition,  even  with  the  less  than 
adequate  ANTC-117  parameters,  would  result  in  an  excessive 

• number  of  alarms,  many  of  which  would  be  unnecessary.  More- 
over, if  the  correct  modified  tau  parameter  values  in  accord- 
ance with  the  theory  of  Ref.  2 were  used  the  expected  number  of 
alarms  would  be  even  greater. 

A CAS  can  reduce  the  number  of  alarms  It  must  give  by  re- 
ducing its  uncertainty  in  estimating  the  location  and  movement 
of  intruding  aircraft.  One  way  it  can  do  this  is  to  use  more 
information  about  the  position  and  motion  of  each  intruder, 
e.g.,  by  adding  bearing  and  bearing  rate  to  the  data  now  re- 
quired for  the  ANTC-117  threat  determination  logic.  A CAS, 
referred  to  as  LCAS,  recently  proposed  by  George  Litchford, 

provides  bearing  data  along  with  the  usual  range  and  altitude. 

1 

In  principle,  the  new  threat  determination  algorithm  deri- 
ved here  makes  the  most  efficient  possible  use  of  the  range. 


By  "efficient,"  as  the  term  is  used  here,  it  is  meant  that  an 
alarm  is  given  if  and  only  if,  in  view  of  the  data  available 
to  the  CAS,  a collision  Is  possible  unless  an  avoidance  maneu- 
ver begins  before  the  advent  of  the  next  measurement  cycle, 
which  is  sometimes  referred  to  as  an  "epoch." 
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altitude  and  bearing  data  assumed  to  be  available.*  Hopefully, 
this  will  result  in  significantly  fewer  unnecessary  alarms  when 
implemented  by  a CAS  than  would  be  the  case  if  the  modified  tau 
condition  were  used  instead. 

The  data  provided  by  LCAS  are  sufficient  to  locate  and 
track  an  intruding  aircraft  with  a precision  that  is  limited 
only  by  measurement  errors.  However,  the  intruder,  whose  inten- 
tions are  assumed  to  be  unknown,  may  deviate  at  any  time  from 
linear  flight,  accelerating  in  the  process,  usually  by  turning, 
as  much  as  1/2  g.  The  threat  determination  algorithm  must 
recognize  a hazardous  situation,  allowing  for  measurement  errors 
and  the  possibility  of  such  an  unknown  acceleration,  in  suffi- 
cient time  to  permit  the  successful  completion  of  an  avoidance 
maneuver . 

Thus,  even  when  bearing  data  are  available  and  the  most 
efficient  possible  threat  analysis  is  used  there  will  still 
exist  enough  uncertainty  concerning  the  intruder's  flight  path 
to  cause  an  alarm  on  some  occasions  when  a collision  would  not 
have  occurred  had  the  alarm  been  ignored.  No  attempt  has  been 
made  here  to  estimate  how  often  this  will  happen  nor  whether, 
in  fact,  the  lower  alarm  rate  generated  by  the  new  algorithm 
will  not  still  be  too  high  for  satisfactory  CAS  performance  in 
heavy  traffic. 

Section  II  discusses  a mathematical  condition  for  the  oc- 
currence of  a mid-air  collision  at  a given  instant  of  tine  when 
the  relative  acceleration  between  the  encountering  aircraft  is 
bounded  by  a known  constant.  Section  III  establishes  an  effi- 
cient alarm  criterion  based  on  that  condition  and  an  assumed 

•This  is  true  if  no  measurement  errors  are  assumed  to  exist. 

When  measurement  errors  are  taken  into  account  in  Section  V 
It  Is  expedient  to  sacrifice  some  efficiency  in  order  to  avoid 
mathematical  complications;  i.e.,  some  alarms  will  be  given 
which  theoretically  might  be  avoided  by  more  effective  use  of 
the  available  data. 


minimum  escape  time,  i.e.,  the  time  required  for  a CAS  to  gener- 
ate the  alarm  and  the  aircraft  to  maneuver  to  safety. 

A threat  algorithm,  which  does  not  take  into  account 
measurement  errors,  will  be  derived  (with  the  aid  of  the  appen- 
dices) and  presented  in  Section  IV.  In  Section  V,  a revised 
threat  algorithm  will  be  given  that  does  include  the  effect  of 
measurement  errors. 


II.  CONDITION  FOR  A MID-AIR  COLLISION 

In  CAS  threat  analyses  (e.g.,  Ref.  2)  it  is  customary  to 
base  the  concept  of  a hazard  upon  the  possibility  of  a collision 
between  two  aircraft  assuming  that  both  obey  regulations  which 
limit  an  aircraft's  speed  and  acceleration.  The  hazard  may 
exist  if  the  collision  can  occur  for  an  allowable  relative 
acceleration  although,  since  it  is  assumed  that  at  least  one 
pilot's  intentions  are  unknown,  the  collision  is  not  actually 
predicted . 

A CAS  officially  recognizes  a hazardous  situation  requir- 
ing an  alarm  when  it  determines  that  after  its  next  observation 
such  a collision  could  occur  before  the  maximum  estimated  time 
(at  some  confidence  level)  needed  for  maneuvering  to  safety  has 
elapsed.  In  this  context  "maneuvering  to  safety"  means  adopt- 
ing a flight  path  (usually  in  a vertical  plane)  which  guarantees 
that  the  encountering  aircraft  remain  separated  by  at  least 
150  ft.* 

• 

Following  Ref.  2,  the  condition  for  a collision  can  be 
stated  mathematically  in  terms  of  the  initial  position  vector 
R of  the  intruder  aircraft  relative  to  the  CAS,  the  initial 
relative  velocity  vector  R,  and  the  relative  acceleration  vector 
A(t).  The  acceleration  A(t)  can  vary  as  a function  of  time  t 
but,  by  assumption,  remains  bounded  in  magnitude.  The  bound  is 
a predetermined  constant  U;  i.e., 

|A(t)|*U 


* 

Cf.  Ref.  1. 
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for  all  t.  The  condition  for  collision  before  escape  is  pos- 
sible is  then 

• f r 

R+tR+l  I A (x)  dids  = o (1) 

~ J O J 0 ~ 

for  some  t in  the  interval 

0 s t < t , 
e 

where  t is  the  estimated  time  required  for  escape.* 

In  Ref.  2 it  was  shown  that  whenever  there  exists  an 
acceleration  A(t)  bounded  in  magnitude  by  a constant  U,  i.e., 

I A ( t )|  £ U, 

for  which  the  collision  condition  (1)  holds,  then  there  exists 
a constant  vector  A also  bounded  in  magnitude  by  U,  i.e., 

s|S  u, 

such  that  (1)  holds  for  some  time  t'  in  the  interval 

0 s t*  * t 

e 

with£(t)  replaced  by  A;  i.e.,  (1)  takes  the  form 

R + t'  R + l/2t'2  A = o.  (2) 

Therefore,  without  loss  of  generality  (2),  which  Is  more  expli- 
cit than  (1),  can  be  used  as  the  collision  condition  for  CAS 
threat  analyses. 


*Not  only  the  time  required  for  any  necessary  avoidance  maneu- 
vers but  also  all  delays  such  as  those  due  to  pilot  and  air- 
craft reactions,  as  well  as  an  allowance  of  one  full  epoch 
for  updating,  are  included  in  tg. 
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Ill . 


THE  CAS  ALARM  CRITERION 


In  order  to  specify  a CAS  alarm  criterion  suitable  for 
translation  into  an  algorithm,  following  Ref.  2 an  attempt  can 
be  made  to  replace  the  vector  collision  condition  (2)  with  a 
scalar  condition  which  does  not  contain  the  unknown  equivalent 
acceleration  vector  £.  The  new  condition  is* 

|R  + te  R|2  < 1/4  t*  U?  (3) 

Indeed,  it  can  be  shown  that  if  (3)  holds  there  exists  a 
t in  the  interval 


0 s t < t 

e 

and  a constant  vector  A,  such  that 

o. 

for  which  the  collision  condition  (2)  is  satisfied,  i.e., 
which 

- • l2 

R + t R + ^ A = o . 

«N»  ^ C ^ 


for 

(A) 


The  proof  is 


not  difficult 
^U2 


b (t ) = 


U — 


Let  b(t)  be  defined  by 

|2 


- - R + t R 


(5) 


*As  Ref.  A has  observed,  a condition  such  as  (3)  is  technically 
not  the  most  efficient  one  possible,  since  it  does  not  include 
the  effect  of  limiting  each  aircraft's  speed,  although  it  may 
be  presumed  that  a maximum  speed  will  be  as  well  enforced  as 
the  maximum  acceleration  which  is  taken  into  account.  As 
Ref.  A has  indicated,  however,  the  inclusion  of  a speed  limit 
should  not  greatly  reduce  the  number  of  expected  alarms  gener- 
ated by  an  otherwise  efficient  threat  algorithm. 


7 


r 


Then  (3)  is  equivalent  to 

b(tg)  > o,  (6) 

which,  because  b(t)  is  continuous,  implies  that  there  is  a t, 
with 

0 s E < t , 
e 

such  that 

b(t ) > 0. 

Then  if  A is  given  by 

A = _ 2 (r  + ^ p), 

it  follows  from  (5)  that 

j&I"  » 

and,  by  substitution,  that  (4)  is  satisfied  as  stated.  Thus, 
if,  for  a given  U and  t , for  some  pair  of  vectors  R and  R (3) 
is  satisfied,  then  a hazard  as  defined  by  (2)  exists.  However, 
the  converse  is  not  necessarily  true.  That  is,  if  condition 
(2)  holds,  so  that  a hazard  exists,  it  does  not  necessarily 
follow  that  (3)  holds.  Therefore,  if  (3)  were  used  alone  as  an 
algorithm  to  determine  the  existence  of  a threat,  in  some  cases 
it  would  fail  to  do  so. 

In  order  to  simplify  the  condition  and  remedy  this  defect 
Ref.  2 replaced  (3)  with  the  modified  tau  criterion 

R + R te  < 1/2  U tg  . (7) 

This  condition  involves  the  scalar  range  R,  which  is  the  magni- 
tude of  the  relative  position  vector  R,  and  the  scalar  range 
rate  R,  which  is  the  component  of  the  relative  velocity  vector 
£ along  £. 

The  modified  tau  condition  (7)  holds  whenever  (2)  holds 
and  therefore  always  predicts  a hazard  when  one  exists.  In 
addition.  Ref.  2 showed  that  when  (7)  is  violated  a hazard  as 
defined  by  (2)  does  not  exi-t.  A proof  of  this  fact  is  also 
sketched  here  in  Appendix  B. 
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However,  a hazard  does  not  necessarily  exist  when  (7)  is 
satisfied.  For  this  reason  that  condition  is  expected  to  pre- 
dict too  many  false  threats  in  heavy  traffic  (cf.  Ref.  3), 
although  Ref.  2 showed  that  (7)  is  the  most  efficient  threat 
determination  criterion  possible,  given  only  a knowledge  of 
R and  R. 

If  the  complete  vectors  R and  £ are  known  it  should  be 
possible  to  replace  (7)  with  a more  efficient  algorithm.  A 
candidate  for  this  role  is  (3);  however,  it  has  already  been 
observed  that  (3)  by  itself  is  not  adequate  because  it  fails 
to  detect  a threat  in  some  cases.  This  defect  can  be  remedied 
by  introducing  some  conditions  to  be  used  in  addition  to  (3). 

It  is  then  possible  to  formulate  a satisfactory  algorithm  which 
detects  all  true  and  no  false  threats  as  they  have  been  defined 
here . 


To  this  end  it  is  useful  to  consider,  first,  a class  of 
cases  for  which  a violation  of  (3)  does  imply  that  no  hazard 
exists.  Ref.  2 proved  that  this  is  always  true  for  the  modified 
t criterion  (7)  (cf.  also  Appendix  B here),  and  that  proof 
depends  on  the  fact  that  the  equation 

l/2Ut2  - Rt  - R = 0 


always  has  exactly  one  positive  root.  If  the  parameters  which 
determine  the  function  b(t)  defined  by  (5)  are  such  that  the 
equation 

b(t)  = 0 (8) 


has  exactly  one  positive  root,  the  same  conclusion  would  follow 
for  the  criterion  (3)-  That  is,  if  a collision  could  occur  at 
a time  t in  the  interval 


i.e . , 


0 s E < t , 
e 

b(E)  - 0, 
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then,  because  b(t)  becomes  arbitrarily  large  as  t becomes  large, 
if  b(t)  cannot  become  zero  again  for  t positive  it  follows  that 

b(te)  > 0 

as  required. 

In  summary,  if  (3)  is  satisfied,  then  an  alarm  should  be 
given,  and  if  (7)  is  violated,  i.e.,  if 

R + Rte  * 1/2  Utg  , (9) 

no  alarm  should  be  given.  It  remains  to  investigate  what 
happens  when  (3)  is  violated  and  (7)  Is  satisfied,  and  the 
answer  depends  upon  the  positive  roots  of  equation  (8). 

In  Appendix  B,  a condition  which  determines  whether  (8) 
has  exactly  one  positive  root  is  derived.  In  the  alternative 
case,  in  which  it  has  more  than  one  positive  root,  special  con- 
ditions to  determine  whether  one  of  those  roots  does,  in  fact, 
lie  between  0 and  t are  also  derived.  Together  with  (3)  and 
(9)  these  additional  conditions  comprise  an  algorithm  which 
satisfies  the  requirements  of  a CAS  alarm  criterion. 

Given  a knowledge  of  the  vectors  R and  R and  the  fact  that 
the  relative  acceleration  between  the  encountering  aircraft  Is 
bounded  by  a specified  constant  U,  but  ignoring  any  speed  limit 
and  measurement  errors,  this  algorithm  provides  the  most  effi- 
cient alarm  criterion  possible.  In  the  next  section  it  is 
given  explicitly  in  terms  of  the  slant  range  R from  the  CAS  to 
the  intruder,  the  range  rate  R,  the  bearing  rate  , the  relative 
altitude  z,  and  the  altitude  rate  z,  all  of  which  are  assumed 
to  be  quantities  known  through  measurement. 
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IV.  THE  THREAT  ALGORITHM 


The  following  parameters  are  assumed  to  be  known  through 
measurement : 

R,  slant  range; 

R,  range  rate; 

4> , bearing; 

$,  bearing  rate; 

z,  relative  altitude  (altitude  difference); 

z,  altitude  range 

These  quantities  are  also  defined  in  Appendix  A in  terms  of 
the  coordinates  of  an  intruder  aircraft  relative  to  a CAS  locat- 
ed at  the  origin  of  the  appropriate  coordinate  system.  In  addi- 
tion, it  is  assumed  that  there  is  a known  bound  U on  the  magni- 
tude of  the  relative  acceleration  between  encountering  aircraft 
and  a known  time  tg,  the  escape  time,  required  for  an  aircraft 
to  maneuver  to  safety  after  it  receives  an  alarm. 


A polynomial  f(t)  in  the  time  variable  t with  coefficients 
determined  by  the  measured  parameters  and  the  quantities  U and 
t is  defined  by  setting 

u2 

b(t)  = J f(t), 

so  that,  in  accordance  with  (5),  f(t)  has  the  form 

4 2 

f(t)  * t - a2t  + a1  t - aQ,  a2io,  a0>o. 

According  to  Appendix  A, 


a 


o 


4R 

U 


(10) 
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a. 


~2  ft 

u L 


Fi  + <T  (R*  - z ; + 


(ftz-zR)^ 

r2  .2 

H - Z 


In  the  absence  of  measurement  errors  the  algorithm  which 
determines  whether  a hazard  exists  and  therefore  an  alarm 
should  be  given  is  derived  in  Appendix  B.  The  algorithm  is 
defined  by  the  following  sequence  of  tests: 


if  yes  then  (2) 


if  no  then  (3)  • 


if  yes  then  (4) 


(1)  R + R t < 1/2  Ut*  ? 

C C 

If  no  then  no  alarm; 


(2)  t*  - a2  t l*  te  - ao  > „? 


If  yes  then  alarm; 


(3)  ax  > o? 

If  no  then  no  alarm; 


(4)  8a23  > 27a12? 


if  yes  then  (5 ) • 


If  no  then  no  alarm; 


(5>  £ h a22? 


If  no  then  no  alarm; 


•It  will  be  observed  that  when  the  intruder  is  directly  over- 
head or  directly  below  the  CAS,  i.e.,  when  R * |z|,  the  param- 
eter a?  is  either  ambiguous  or  infinite.  It  may  also  be  ex- 
pected that  a?  will  be  very  sensitive  to  error  when  R~|zj. 
Therefore,  as  a practical  matter  unless  R differs  from  |zj  by 
more  than  the  maximum  combined  error  of  the  measured  parameters 
R and  z the  usual  modified  tau  criterion  should  be  used  instead 
of  the  threat  algorithm.  This  is,  in  fact,  equivalent  to  let- 
ting a2  become  infinite  in  the  threat  algorithm. 
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if  yes  then  (6 ) . 


(6)  for  sQ  = 


12ao  M 
a 2 ' J’ 


al  * 
n72 
2 


2so 


7 


If  no  then  no  alarm; 

if  yes  then  (7 ) . 

(7)  6tei 2  > a2? 

If  yes  then  alarm; 

if  no*  then  (8). 

(8)  4tg3  - 2a2  te  + a1  < 0?*« 

If  yes  then  alarm; 

if  no  then  no  alarm. 

Example  1 . 

As  an  example  of  how  this  algorithm  would  be  applied,  sup- 
pose that  the  following  CAS  parameters  are  measured  for  an 
intruder : 


range  R * 4000  ft, 
range  rate  R = -590  ft/sec, 
relative  altitude  z ■ 200  ft 
relative  altitude  rate  z ■ -1600  ft/min  = -26.67  ft/sec 
bearing  rate  <j>  = l°/sec  ■ 0.01745  rad/sec 
If  the  maximum  relative  acceleration  is  assumed  to  be  lg,  then 

U = 32.2  ft/sec2. 


i ^ 

6t  / a2  because  of  steps  (4)  and  (5). 

• * 3 

4tg  - 2a2tg  + a^  ^ 0 because  of  steps  (6)  and  (7). 
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Then  from  (10)  it  is  found  that 


aQ  * 6.1726  x 10*4  sec14, 
ax  * 1.8209  x lO4  sec3, 

p 

a2  = 1361.7  sec  . 


If  it  is  assumed  that  the  escape  time  t is  25  sec,  then 
at  step  (1)  of  the  algorithm  it  is  found  that 

R + Rt  = -10750  ft  < 10,062  = H»Ute2. 

Thus,  an  alarm  would  be  given  by  the  modified  tau  criterion, 
but  for  the  threat  algorithm  of  this  section  it  would  be  neces- 
sary at  this  point  to  consider  step  (2).  Then  it  would  oe 
found  that 

t 14  - a«t + a.t  - a.  = -6.6938  x 1014  sec14  < 0. 
e 2 e 1 e o 


Thus,  it  would  be  necessary  to  consider  step  (3): 

a^  = 1.8209  x 10*4  sec3  > o, 

which  leads  to  step  (iJ).  It  is  found  that 

8a23  * 2.0199  x 1010  > 8.9523  x 109  * 27a12, 
which  leads  to  step  (5).  For  that  test 


a = 6.1726  x lO*4  < 1.5^52  x 105  = ^ a 2 

O 1 .C.  d. 

Thus,  it  is  necessary  to  go  to  step  (6).  For  this  purpose  s 
must  be  calculated: 


iAF fWi 


0.19367 


Then  step  (6)  gives 

a. 


a. 


172 


0.36238  i 0.35828 


2s  - 1*3  - 
o o 
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Thus,  according  to  the  test  at  step  (6)  step  (7)  must  be  con- 
sidered. At  step  (7)  it  is  found  that 

6te2  = 3750  > 1361.7  = a2; 
hence  an  alarm  must  be  given. 

This  result  can  be  verified  by  inspection  of  Fig.  1 which 
shows  the  curve 

f(t)  = - a2t2  + a1t  - aQ. 

It  may  be  observed  that  f(t)  vanishes  at  t = 6 sec  and  again 
at  a little  more  than  8 sec,  both  of  which  represent  possible 
collision  times  before  tg  = 25  sec. 

Example  2. 

The  purpose  of  this  example  is  to  see  what  happens  if  the 
bearing  rate  is  changed  but  all  of  the  other  parameters  of  ex- 
ample 1 are  retained.  The  modified  tau  criterion  is  unaffected 
by  such  a change.  However,  if  the  bearing  rate  4>  is  taken  to 
be  2°/sec  (0.034907  rad/sec),  then  aQ  and  a^  remain  the  same 
but  now 

a2  = 1418  sec2. 

In  this  case  it  will  be  found  that  the  threat  algorithm  indi- 
cates "no  alarm"  at  step  (6). 

Indeed,  Fig.  1,  which  contains  the  curve  for  f(t)  shows 
that  f(t)  has  a maximum  at  about  7 sec,  is  negative  between 
0 and  25  sec,  and  has  positive  curvature  at  25  sec.  Appendix 
B demonstrates  that  these  facts  should  account  for  the  indi- 
cated behavior  of  the  threat  algorithm.  Since,  as  the  curve 
shows,  f(t)  does  not  vanish  between  0 and  25  sec  there  is  no 
threat  of  collision  before  the  escape  time  has  elapsed. 

Example  3. 

If  instead  of  the  bearing  rate  4>  changing,  it  is  the  range 
which  is  changed  to  3 nmi  (18,240  ft)  while  the  other  parameters 
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remain  the  same  (e.g.,  $ stays  l°/sec),  the  threat  algorithm 
indicates  an  alarm  at  step  (2).  That  is,  the  coefficients 
obtained  from  (10)  are 

aQ  = 1.2835  x 108  sec14, 

ax  = 8.3034  x 104  sec3, 

&2  = 1735-4  sec2, 

so  that 

t e ) = te4  - a^te2  + a1te  - aQ  - 9.8320  x 10*4  > 0. 

Fxample  4. 

On  the  other  hand,  if  the  range  R is  increased  to  4 nmi 
(24,320  ft)  with  the  other  parameters  left  unchanged,  then 

aQ  = 2.2818  x 108  sec14, 

= 1.1071  x 105  sec3, 

= 2039-8  sec2, 

and  at  step  (4)  it  is  found  that 

8a23  = 6.7897  x 1010  < 3.3093  x 1011  = 27a12. 

Thus,  no  alarm  would  be  generated  for  the  larger  range.  Fig.  2 
shows  the  curves  for  f(t)  corresponding  to  the  3 and  4 nmi 
ranges.  Neither  has  a maximum  (as  expected,  since  both  fail 
the  test  at  step  (4)),  the  3 nmi  curve  is  zero  at  a time  slightly 
beyond  23  sec,  indicating  a collision  threat  before  the  escape 
time  of  25  sec,  and  the  4 nmi  curve  remains  negative  between  0 
and  25  sec,  indicating  no  threat. 
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V.  THREAT  ALGORITHM  WITH  MEASUREMENT  ERRORS 

A.  ERROR  TOLERANCES  AND  PARAMETER  UPPER  AND  LOWER  BOUNDS 

Each  of  the  parameters  measured  by  a CAS  will  have  some 
specified  error  tolerance,  so  that  a single  measurement  will 
actually  define  a range  of  possible  values  for  a parameter 
rather  than  a single  value.  Therefore,  it  is  necessary  to  re- 
place the  algorithm  given  in  Section  IV  by  one  which  will 
generate  an  alarm  whenever  a measurement  and  the  error  toler- 
ances are  consistent  with  any  set  of  independent  parameter 
values  that  satisfy  an  alarm  condition  of  the  original  algo- 
rithm. 

If  the  revised  algorithm  were,  indeed,  restricted  to 
indicating  an  alarm  only  in  such  a circumstance  its  efficiency 
would  remain  unimpaired,  although  the  number  of  alarms  would 
necessarily  increase  because  of  the  contribution  of  measurement 
error  to  the  uncertainty  in  estimating  an  intruder's  position 
and  movement.  However,  it  will  be  found  mathematically  con- 
venient to  widen  the  effective  tolerance  ranges  of  certain 
derived  quantities,  replacing  their  possible  maximum  and  mini- 
mum values  with  grosser  estimates  which  are  no  more  than  upper 
and  lower  bounds.  This  will  reduce  the  algorithm's  efficiency 
somewhat,  thereby  Increasing  the  number  of  alarms  it  will  gene- 
rate. The  result  is  equivalent  to  an  increase  in  the  error 
sensitivity  of  the  affected  quantities. 

In  revising  the  threat  algorithm  of  Section  IV  it  will  be 
convenient  to  use  the  following  notation.  When  the  subscript 
m is  added  to  the  symbol  for  a measured  quantity  the  new  sym- 
bol will  designate  the  quantity's  minimum  value  as  determined 
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by  the  specified  error  tolerance.  Similarly,  adding  the  sub- 
script M will  provide  a symbol  for  the  quantity's  maximum. 
Thus,  for  example,  if  a range  measurement  is  R and  the  error 
tolerance  is  ± e,  so  that  the  range  R would  be  properly  in- 
dicated by 


then 


and 


R = H0  ± e. 


R = R - e 
m o 


rM  = Ro  + e * 


where  the  usual  convention  that  e is  non  negative  has  been 
assumed .* 

On  the  other  hand,  for  symbols  representing  quantities 
which  are  derived  rather  than  measured  directly,  the  addition 
of  the  subscripts  m and  M will  indicate,  less  precisely,  just 
lower  and  upper  bounds.  These  may  or  may  not  be  the  true  mini- 
ma and  maxima  implied  by  the  given  measurement  error  tolerances. 

For  example,  the  quantity  aQ  defined  by  (10)  will  lead 
naturally  to  the  definitions 


4r_ 


a 


m 


and 


om 


o M 


U 


4R 


M 


since  R is  non  negative  by  virtue  of  its  geometrical  meaning. 

However,  for  a^,  defined  by  (10),  it  would  be  permissible,  for 

*TEF-  error  e Is  a fixed  bound  on  the  deviation  of  R from  R0. 

A statistical  error  fluctuation  would  be  taken  into  account 
by  assigning  some  kind  of  confidence  limit  along  with  e,  e.g., 
by  stating  that  e bounds  the  error  with  99%  probability  or 
that  e is  a 3o  bound  on  the  error. 
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example,  (although  not  recommended  in  this  case)  to  define 


aim  = ~8  RN! 


M 


which  is  a lower  bound  smaller  than  the  minimum  of  a-,  whenever 

V 0.* 

A known  physical  or  geometrical  relationship  may  play  a 

role  in  determining  the  minimum  or  maximum  value  of  a derived 

2 

quantity’s  tolerance  interval.  For  example,  the  quantity  p 
defined  in  terms  of  the  relative  range  R and  relative  altitude 
z by 

2 2 2 

p = R^  - z 

must  always  be  positive,  whatever  the  tolerance  limits  may  be 

? 

on  R and  z separately.  Thus,  since  p cannot  be  less  than  zero 
it  follows  that 


2 2 , i 2 *» 

p = R - 
m m 


‘M 


(ID 


if 


m 


Z M * ° 


but 


p = o 
m 


(12) 


if 


m 


Z M s °* 


*In  this  case  the  true  minimum  value  of  a1  would  be  -8RmRM  > 


which  is  positive. 


U 


«* 


Since  z can  be  either  positive  or  negative  |z 
to  z^2  if  Zjyj  is  positive;  if  zM  is  negative, 


M 


is  only  equal 


M 


Jm 
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On  the  other  hand,  it  is  always  understood  that 

4 - V1  - !*!,/■  <»> 

B.  AUXILIARY  SUBALGORITHMS 

Before  attempting  to  revise  the  threat  algorithm  of  Section 
IV  it  will  be  found  useful  to  formulate  some  auxiliary  algo- 
rithms which  will  help  to  simplify  the  description  of  the 
necessary  logic.  These  cover  the  arithmetical  steps  which  are 
implied  when  the  subscripts  m and  M are  added  to  some  simple 
algebraic  expressions. 

The  addition  of  a subscript  to  an  expression  of  the  form 
(u  v)  or  |u  + v|  has  the  effect  of  operating  on  the  expression 
to  produce  an  algebraic  function  of  the  quantities  u , u^  vffll 
and  Vj^  Four  such  operations  are  possible:  (u  v)m,  (a  v)^, 

|u  + v|m,  and  |u  + vl^.  The  results,  which  may  be  regarded  as 
algorithms  for  the  calculation  of  the  subscripted  quantities, 
are  given  in  Tables  1,  2 and  3. 

It  should  be  emphasized  here  that  these  results  can  be 
relied  upon  to  provide  actual  maxima  and  minima,  in  general, 
only  when  u and  v are  independent.  If  u and  v are  dependent 
inefficient  bounds  may  be  given  instead,  e.g.,  in  the  case 

p 

(u  ) which  can  be  calculated  from  Table  2 by  setting  u = v. 
m ^ 

When  u < 0 < uM,  Table  2 gives  a lower  bound  for  (u‘‘)  below 
m P'l  2 2 ^ 

the  true  minimum,  which  is  min  (u  , u^ ) . 


TABLE  1.  Subalgorithm  for  I u + v | m and  |u+v|M 


if  < 0 

If  U+Vm 

m m 

> 0 

If  U+V£0*Um+Vm 
mm  mm 

lu  + VL 

1 m 

-<UM  * VM> 

um  + 

vm 

0 

lu  + VIM 

-(um  * v»> 

> + 

VM 

max  [•<VvJ-(uM+vM>] 

TABLE  2. 


Subalgorithm  for  (uv)m 


In  addition  to  the  algorithms  defined  by  the  tables, 
formulas  for  the  expressions  (^)  and  (^)M>  where  u and  v are 
both  non  negative,  are  needed.  These  are 


and 


(I'D 


(15) 


Most  of  the  subscripting  required  for  the  revised  threat 
algorithm  can  be  accomplished  by  repeated  applications  of  the 
subalgorithms  defined  by  Tables  1,  2,  and  3,  the  formulas  (14) 
and  (15),  and  the  trivial  relations 


(u  + v ) “ u + v . 

m m m 


(u  + V)M  * UM  * VM’ 
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It  should  be  observed,  however,  that  this  procedure  is  not 
uniquely  defined  for  a given  expression  and,  therefore,  may  lead 
to  different  results,  depending  upon  the  grouping  of  terms  and 
the  order  in  which  the  subscripting  is  done. 

If  at  each  application  of  a subalgorithm  in  this  procedure 
u and  v are  independent,  the  true  minimum  or  maximum  of  the 
original  expression  will  be  obtained.  If,  however,  at  some 
step  u and  v are  dependent  the  final  result  may  only  be  a bound, 
i.e.,  a quantity  which  is  smaller  than  the  true  minimum  or 
larger  than  the  true  maximum.  As  an  example,  consider  the  cal- 

p 

culation  of  (x  + x)^,  where 


If  the  calculation  is  done  as  follows: 

(x2  + x)M  = (x2)m  + xM  = 4 + 0 = 4, 

the  result  is  not  a true  maximum  but  an  upper  bound.  In  fact. 
It  is  easy  to  see  that  the  true  maximum  occurs  when  x * x = -2 
and  that  its  value  is  2. 

The  calculation  can  also  be  done  by  first  setting 
(x2  + x)M  = ( x ( x + 1))M, 

next  Identifying 


u = x,  v = x + 1 , 

and  finally  using  Table  3.  It  will  be  found,  first,  that  in 
this  case 


and 


m 


-2, 


M 


0 


m 


-1, 


M 


1. 


Table  3 applied  to  (uv)M  will  then  give  the  true  maximum  2 as 
the  result. 
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The  calculation  of  the  extrema  a , a a.  , a,  a„  , 

om*  oM’  lm’  1M*  2m’ 

a2M  corresponding  to  the  coefficients  given  by  (10)  could  be 
accomplished  by  repeated  applications  of  the  subalgorithms  but 
with  some  unnecessary  loss  of  efficiency.  Instead,  it  is  pre- 
ferable to  use 


a 

om 


aoM 


(16) 


for  a and  a . Tables  2 and  3 can  be  used  to  calculate  a. 
om  oM  lm 

and  a^M  without  loss  of  efficiency,  however. 


For  a2m  and  a2M,  accepting  some  loss  of  efficiency  is 
probably  worthwhile  for  the  sake  of  mathematical  simplicity  in 
the  calculations.  A reasonable  form  of  a2  for  the  calculation 
is 


*2  = i_  [|R|2  + l*l2  p2  + 

U2  L 2 — 

o 


so  that 


a. 


'2m 


U1* 

U L 


Rlm2  + l*lm2  »m2  + H 


PK 


(17) 


and 


2M 


[iRl 

7 L 


M2  + UIm2  PM2  + ^ 


^|RIm2] 

n -J 


m 


(18) 


In  (17)  and  (18),  the  second  and  third  terms  can  be  calculated 

from  Tables  1,  2,  and  3 with  the  aid  of  (11),  (12),  and  (13) 

2 2 

for  calculating  pm  and  pM  . 

Expressions  for  the  extrema  sQm  and  sqM  of  the  quantity 
sQ,  which  was  defined  in  step  (6)  of  the  original  threat  algo- 
rithm of  Section  IV,  can  be  obtained  by  inspection.  For  this 
purpose  it  should  be  recalled  that  aG  and  a2  are  both  non  neg- 
ative. Although  some  loss  of  efficiency  occurs  because  aQ  and 
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a2  are  dependent,  reasonable  expressions  are 


f 


I 


I 


I 


t 


I 


1 


and 


(19) 


(20) 


It  will  be  seen  later  that  the  values  of  s and  s given  by 

om  oM 

(19)  and  (20)  are  guaranteed  to  be  real  because  the  revised 
threat  algorithm  would  otherwise  have  terminated  at  step  (5) 
before  a need  for  these  quantities  could  have  arisen. 


Algorithms  are  also  needed  for  the  extrema  f(t  ) , f(t  )M, 

0 III  6 n 

f'(t^)™>  and  f' (t  )M  of  the  polynomials  f(t  ) and  f'(t  ) which 

0 m 0 M 0 0 

appear  in  steps  (2)  and  (8)  of  the  threat  algorithm  of  Section 
IV.  The  original  polynomials  were  defined  by 


f(te) 


a t 
2ce 


a.  t 

1 e 


a 

o 


(21) 


and 


f(te)  - - 2a2te  * ar 

For  f'(t  ) it  is  probably  advisable  to  accept 
0 

f ’ (t  ) = 4t  3 - 2aOMt  + a, 

e m e 2M  e lm 


and 


f'(t  )M  * 4t  3 - 2a0mt  + a,M. 
e M e 2m  e 1M 


(22) 


(23) 

(24) 


However,  somewhat  more  efficient  bounds  can  be  obtained  for 
f(t&)  by  first  expressing  it  in  the  form,  using  (10), 

f(t  ) - o - \ |R  + t R|2,  (25) 

e u 


where 


2 , Rz-zR 
) + 

P 


(26) 


w 
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wherein  further  subscripting  can  be  carried  out  with  the  aid  of 
the  previously  derived  subalgorithms.  The  subalgorithms  can  be 
applied  along  with  (27)  and  (28)  to  (25)  in  order  to  calculate 
f(te>m  and  f(te)„. 

C.  THE  REVISED  THREAT  ALGORITHM 

In  order  to  obtain  the  revised  threat  algorithm  it  is  only 
necessary  to  strengthen  the  Inequality  condition  for  no  alarm 
or  weaken  the  condition  for  an  alarm  at  each  step  of  the  threat 
algorithm  in  Section  III  by  applying  the  appropriate  extremum 
operations.  The  subalgorithms  derived  in  Section  VB  now  make 
this  a completely  mechanical  procedure. 

The  result  is: 

(1)  R + R t < i U t 2? 

m me  d e 

If  no  then  no  alarm; 

if  yes  then  (2) . 

(2)  f(te)M  > 0? 

If  yes  then  alarm; 

if  no  than  ( 3) • 

(3)  a1M  > 0? 

If  no  then  no  alarm; 

if  yes  then  (4 ) . 

(4)  8a2M  > 27alm? 


I 
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if  yes  then  (5) . 


if  yes  then  (6) . 


if  yes  then  (7) 


if  no  then  (8) 


If  no  then  no  alarm; 


(5)  aom  ‘ h Vn  ? 

If  no  then  no  alarm; 


( 6 ) a,  M 2 2s  - 4 s „ ? 
1M  om  oM 

1372 

m 

If  no  then  no  alarm; 


(7>  6te  > a2m! 


If  yes  then  alarm: 


(8)  f'(te)m  < 0? 


If  yes  then  alarm; 


if  no  then  no  alarm. 


For  an  example  to  illustrate  the  use  of  the  revised  threat 
algorithm,  the  following  error  tolerances  will  be  assumed: 


range 

(R) 

± 

200  ft. 

range  rate 

(R) 

± 

50  ft/sec  (about  ± 30kn) 

altitude 

( z) 

t 

200  ft. 

altitude  rate 

(z) 

± 

56.6  ft/sec 

bearing  rate 

<♦) 

t 

5 /sec. 

Then,  for  measured  parameters  having  the  values 

R - 12160  ft  (2  nmi ) , 

A ■ -590  ft/sec, 
z - 200  ft, 
z - -26 . 67  ft/sec , 
i ■ l6°/sec, 
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» 


a 


j 


I 


the  extrema  which  must  be  used  in  the  threat  algorithm  will 
be  given  by 

Rm  = 11,960  ft,  Rm  = 12,360  ft, 

Rm  = -640  ft/sec,  Rm  = -540  ft/sec, 


z = 0 ft, 
m 


zM  * 400  ft. 


zm  = -83.3  ft/sec,  zM  = 29.9  ft/sec. 


11  /sec , 


<t>M  = 21  /sec 


Using  the  methods  described  in  this  section  for  calculat- 
ing the  derived  parameter  extrema,  it  will  be  found  that 

Pm2  = 1 . 4288x10®  ft,2  pM2  = 1.3277xl08  ft,2 

aom  = 5-5l84xl0^  sec4  aQM  = 5-8937xl03  sec4, 

alm  = 4 • 9831x10^  sec?  a1M  = 6.1035xl04  sec3, 

a2m  = 2.l442xl04  sec2  a2M  = 8.0782xl04  sec2, 

sQm  = 9 • 1970x10” 3 , sQM  = 3.5873x10"2 

f(t  )M  = 1.2313xl07sec4,  f* (t  ) = -3.9268xl06sec3. 

e M * em 

The  threat  algorithm  applied  to  these  quantities  then  leads  to 
the  conclusion  at  step  (8),  that  an  alarm  should  be  generated. 


that 


If  the  measured  bearing  rate  $ is  increased  to  17  /sec,  so 


i = 12°/sec  and  = 22°/sec, 
m m 


while  the  other  measured  parameters  are  left  unchanged,  then 

2 2 

the  derived  parameters  Pm  > PM  . aom»  a0M’  alm’  and  alM  remain 
unchanged  but  new  values  for  a2m,  a2M  som>  and  soM  are  obtained, 
These  are  given  by 

a2m  = 2 . 5304xl04sec2 , a2M  * 8 . 8502xl04sec‘i , 
s __  = 8 . 3946xl0~3 , s.w  « 3.038lxl0T2 
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On  applying  the  threat  algorithm  to  the  new  values  of  the  de- 
rived parameter  extrema  it  will  be  found  that  "no  alarm"  is  in- 
dicated at  step  (6). 
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APPENDIX  A 


COEFFICIENTS  IN  THE  THREAT  ALGORITHM 

The  condition  (3)  of  Section  II  is  the  basis  for  the 
hazard  alarm  algorithm  given  in  Section  IV.  That  algorithm 
operates  on  the  coefficients  aQ,  a^,  and  a 2 of  a polynomial 
f(t)  derived  from  (3),  which,  in  turn,  depends  upon  the  com- 
ponents of  the  range  vector  R between  the  CAS  and  an  intruder 
aircraft  and  upon  the  time  derivative  R of  R.  The  components 
of  these  vectors  depend  upon  a set  of  quantities  measured  by 
the  CAS: 

R,  the  range  between  the  CAS  and  the  intruder, 

R,  the  range  rate, 

$ , the  bearing  of  the  intruder 
<ji,  the  bearing  rate, 

z,  the  difference  in  altitude  between  the  intruder 
and  the  CAS, 

z,  the  relative  altitude  rate. 

Expressions  for  the  coefficients  a^,  a^ , and  a^  in  terms  of  the 
measured  quantities  will  be  derived  in  this  appendix. 

It  will  be  convenient  to  use  a right-handed  spherical 
coordinate  system  centered  at  the  CAS.  The  polar  axis  is 
perpendicular  to  the  ground.  The  polar  angle  0 is  zero  along 
the  axis  in  the  upward  direction  and  increases  downward  to 
its  extreme  value  of  180°  along  the  axis  toward  the  ground. 

The  azimuth  angle  4>  measures  the  absolute  bearing  of  the  intruder 
and  is  therefore  zero  at  some  predetermined  geographical 
direction,  e.g.,  magnetic  north.  Looking  down,  a horizontal  vector 
which  emanates  from  the  CAS  and  rotates  counterclockwise  will 
have  an  azimuth  <t>  that  Increases  as  it  rotates.  The  radial 
coordinate  R is  the  intruder  range.  A corresponding  vector 
coordinate  system  is  based  upon  a set  of  three  orthogonal  unit 


vectors  aQ , g^,  and  gR  pointing  in  the  direction  of  increasing 
0,  <p , and  R. 

Condition  (3)  of  Section  II  depends  upon  the  vectors  g 
and  g.  The  vector  g is  given  by 

R = R aR  (A-l ) 

in  terms  of  the  spherical  coordinate  system.  The  vector  R is 
given  by 

R = R SR  + R aR.  (A-2 ) 

The  unit  vector  a changes  direction  from  point  to  point  and 

~ n 

therefore,  because  of  aircraft  motion,  it  does  not  remain 
constant  in  time.  For  this  reason,  as  indicated  in  (A-2),  its 
time  derivative  aR  is  not  zero  in  general. 

For  the  purpose  of  calculating  aR  it  is  helpful  to  use 
a Cartesian  coordinate  system,  for  which  the  unit  orthogonal 
basis  vectors  a^ , a^ t and  a^  are  constant  throughout  space 
and,  thus,  do  not  change  in  time.  A right-handed  Cartesian 
coordinate  system  for  this  purpose  may  be  defined  by  letting 
the  unit  vector  a^  point  in  the  direction  0 = o,  the  unit 

vector  a point  in  the  direction  4>  = o in  the  plane  0 = 90°, 

^ o 

and  the  unit  vector  a point  in  the  direction  4>  = 90w  in  the 

o * 

plane  0 = 90  . The  coordinate  systems  are  illustrated  in 
Fig.  A-l. 

The  relationship  between  the  basis  vectors  of  the 
spherical  coordinate  system  and  those  of  the  Cartesian  coordinate 
system  is  given  by  the  following  equations: 


aR  = (cos  ip  sin  6)  gx  + (sin  <p  sin  0)  + (cos  0)  az  , 

= (-sin  4>)ax  + (cos  4> )a.y  , (A-3) 

aQ  = (cos  <p  cos  0)ax  + (sin  <J>  cos  0)a  - (sin  0)a  , 

and 

a - (-sin  4))a^  + (cos  <t>  cos  0)an  + (cos  $ sin  0)aD  , 

X <J>  ~ U ~ n 

( A-4 ) 

a * (cos  4> ) a . + (sin  <p  cos  0)aQ  + (sin  <p  sin  0)§D  , 

-y  ~o  rt 

a = (-sin  0)a„  + (cos  6)  aD  . 

- Z ~ o n 

It  follows  from  (A-3)  and  (A-4)  that 

* / * * 

~R  = s*n  ^ s*n  ® + ® cos  ^ cos  ®)8X 

+ (4>  cos  <{>  sin  0 + 0 sin  0 cos  0)ay 
- (0  sin  0)az  = (<}>  sin  0)g.  + 0^  . 


Then,  with  the  aid  of  this  result,  (A-2)  can  be  written 

R = R a + (R  * sin  0 )a  + Re  a.  (A-5) 

r - <p  u 

• • 

Actually,  while  <p , <p , R,  and  R are  assumed  to  be  measured 
quantities,  6 and  6 are  not  measured  directly.  However,  z 
and  z are,  and  0 and  0 can  be  expressed  in  terms  of  z,  z,  R, 
and  R.  The  pertinent  relationship  is  given  by 

z * R cos  0 , (A-6) 

from  which  it  follows  that 


A-5 


i 


cos  0 


and 


sin  0 


R 


(A-7 ) 


Then,  from  (A6), 


z = ft  cos  0 - R0  sin  0, 


and,  from  ( A 7 ) 

A _ Rz  - Rz 

‘ "ri~ ? ' 

RyR  -z 

If  these  relations  are  substituted  into  (A-5)  the  result  will  be 
an  expression  for  R , 

R = R aR  + (♦  *9  ’ (A'8) 

in  terms  of  measured  quantities  alone. 

The  condition  (3)  of  Section  II  expanded  in  terms  of 
vector  dot  products  is 

uV 

R2  + 2(B  • R)te  + |p|2t2  < -rpS-  . (A-9) 

This  inequality  can  be  written  in  the  form 

f(tg)  > o , 


A- 6 


where 


with 


it  p 

f(t)  = t - a^t  + a^t  - aQ 


a = 


a,  * - 


8(R*ft) 


4 |R 


tr 

2 


a„  = 


U2 


With  the  aid  of  (A-l)  and  (A-8),  remembering  that  the  basis 
vectors  have  unit  magnitude  and  are  mutually  orthogonal,  the 
coefficients  of  f(t)  can  then  be  written  in  terms  of  the  measured 
quantities  alone: 


a2 


\ R2+j>2  (R2-z2)  + 

IT  . 


(Rz-Rz ) 
R -z 


(A- 10) 


C 
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APPENDIX  B 


NECESSARY  CONDITION  FOR  A COLLISION  THREAT 

A.  SUFFICIENT  CONDITIONS  FOR  A THREAT  AND  FOR  NO  THREAT 

Relation  (3)  of  Section  III  is  a sufficient  condition, 
depending  on  range,  altitude,  and  bearing  data,  for  a colli- 
sion threat.  The  condition  is  also  equivalent  to 

f*  (te)  > 0, 

where  f(t)  is  the  biquadratic  polynomial 

f(t)  = t11  - a,  t2  + a,  t - a (B-l) 

c 1 O 

defined  in  Section  IV.  The  coefficients  a , a, , and  a^  are 
given  in  terms  of  the  measured  data  by  (10)  of  Section  IV  or 
( A-10 ) of  Appendix  A.  From  their  definitions  it  is  clear  that 

a„  * o and  a„>  o, 

2 0 

while  can  be  positive,  negative,  or  zero. 

On  the  other  hand,  according  to  Ref.  2,  if  the  modified 
tau  condition  (7)  is  violated,  i.e.,  if 

• (B-2) 
e 

no  threat  exists.  This  can  be  seen  from  the  argument  of 
Section  III  applied  to  a function  b(t)  which  is  redefined  for 
this  purpose  by 

b(t)  = \ Ut2  - Rt  - R. 

The  result  then  follows  from  the  fact,  remarked  In  Section  III, 
that  the  equation 

b ( t ) = j Ut2  - Rt  - R = o 
B-2 


always  has  exactly  one  positive  root.  This  Is  because  the 
quadratic  polynomial  b(t)  is  negative  at  t = o and  increases 
without  limit  as  either  t or  -t  becomes  large.  Since,  accord- 
ing to  (B-2),  b(te)  is  negative  the  positive  root  must  occur 

at  a time  later  than  t . 

e 

The  question  now  left  to  be  answered  in  this  appendix  is 
what  happens  when  neither  (3)  nor  (B-2)  are  satisfied.  That 
is,  what  can  be  concluded  about  the  existence  of  a threat  if 

f(te)  * 0?  ( B—  3 ) 

B.  CRITICAL  POINTS  AND  ROOTS 

It  was  shown  in  Section  III  that  when  the  equation 

f(t)  = o (B-4 ) 

has  only  one  positive  root  and  the  condition  ( B— 3 ) is  satisfied, 

no  collision  can  occur  at  a time  earlier  than  t . Thus,  in 

6 

order  to  determine  whether  an  alarm  should  be  given  when  ( B—  3 ) 

Is  satisfied,  a first  step  might  be  to  find  out  if  (B-4)  has 

more  than  one  positive  root.  If  not,  the  question  will  be 

answered  in  the  negative;  otherwise,  it  will  be  necessary  to 

determine  whether  one  or  more  of  the  roots  occur  for  values  of 

t between  0 and  t . 

e 

Since  f(t)  is  negative  at  t = 0 and  becomes  arbitrarily 
large  as  either  t or  -t  becomes  large,  (B-4)  has  at  least  one 
positive  and  one  negative  root.  It  may,  of  course,  have  as 
many  as  four  distinct  real  roots. 

It  also  follows  that  f(t)  must  have  at  least  one  minimum. 
This  corresponds  to  the  fact  that  the  derivative  f'(t),  being  a 
cubic  polynomial,  must  vanish  for  at  least  one  real  value  of  t. 

If  f'(t)  vanishes  for  only  one  real  value  t of  t then 
the  minimum  of  f(t)  must  occur  at  t and  there  can  be  no  other; 
nor,  for  that  matter,  can  f(t)  have  a maximum  in  such  a case. 
Furthermore,  (B-4)  will  then  have  just  the  one  positive  and  the 

B-3 


one  negative  root.  This  situation  is  illustrated  in  Fig.  3 in 
which  the  minima  of  both  curves  occur  at  negative  values  of  t. 

On  the  other  hand,  if  f'(t)  vanishes  at  three  distinct 
values  of  t,  f(t)  will  have  two  minima  separated  by  a maximum. 
This  relationship  Is  topologically  clear;  It  can  also  be  veri- 
fied by  considering  the  sign  changes  of  the  second  derivative 
f"(t). 

From  ( B— 1 ) , 

f"(t)  = 12t2  - 2a2.  ( B— 5 ) 

Thus  f"(t)  vanishes  at  two  values  t + of  t given  by 

t±  - ±J^_  . (B-6) 

These  are  the  inflection  points  of  f(t)  and  are  evidently 
located  symmetrically  about  zero. 

It  can  be  seen  by  inspection  of  (B-5)  that  the  slope 
f"(t)  of  f'(t),  where  from  (B-l) 

f ' (t)  = 4t3  - 2a2  t + ax,  ( B-7 ) 

is  negative  for  t between  t_  and  t+  and  is  positive  for  all 
other  values.  Moreover,  f'(t),  itself,  has  a maximum  at  t 
and  a minimum  at  t+ . Thus,  if  f'(t)  vanishes  at  three  points, 
one  lies  to  the  left  of  t_  and  must  correspond  to  a minimum  of 
f(t),  one  lies  between  t_  and  t+  and  must  correspond  to  a maxi- 
mum of  f(t),  and  the  third  lies  to  the  right  of  t+  and  must 
correspond  to  another  minimum  of  f(t). 

Since 

f'  (0)  * a2 

it  is  also  clear  that  if  a^  is  negative,  the  maximum  of  f(t) 
lies  between  t and  0,  because  f'(t)  will  change  sign  between 
these  two  values  of  t.  Similarly,  if  a^  is  positive,  the 


maximum  of  f(t)  lies  between  0 and  t,.  This  situation  is 
illustrated  in  Figs.  1 and  2. 

Since  f(0)  is  negative,  it  then  follows  that  for  a^  nega- 
tive, (B-M)  can  only  have  one  positive  root.  For  a1  equal  to 
zero  the  maximum  of  f(t)  will  be  at  zero  and  will  therefore  be 
negative.  Thus,  f(t)  can  have  more  than  one  positive  root  only 
if  a^  is  positive.  It  follows  that  when  ( B— 3 ) is  satisfied  a 
hazard  condition  can  only  exist  when 

a.^  > 0 . 

This  condition  is  equivalent  to  requiring  a negative  value  for 
R,  which  would  imply  that  the  intruder  is  approaching. 

C.  THE  CUBIC  DISCRIMINANT 

Since  f(t)  cannot  vanish  at  more  than  two  points  if  it  has 
just  one  minimum  and  lacks  a maximum,  a necessary  condition 
that  (B-4)  have  more  than  one  positive  root  is  that  f'(t)  vanish 
at  more  than  one  real  value  of  t . In  fact,  the  condition  may 
be  stated  more  restrictively : the  equation 

f'(t)  = 0 (8-8 ) 

must  have  three  distinct  real  roots.  This  follows  from  the 
observation  that  a multiple  root  of  (B-8)  must  occur  at  an 
Inflection  point  of  f(t),  and  at  such  a point  f(t)  will  be 
stationary  but  will  have  neither  a maximum  nor  a minimum. 

The  condition  that  f’(t)  have  three  real  roots  could  be 
obtained  from  the  well  known  cubic  discriminant  (cf.  Ref.  B-l, 
p.  432).  However,  it  may  be  instructive  to  derive  it  here 
directly  in  terms  of  the  coefficients  in  f(t). 

It  has  been  observed  that  f'(t)  has  two  critical  points 
t+,  given  by  (B-6).  It  is  easily  verified  that  t corresponds 
to  a maximum  of  f(t)  and  t to  a minimum.  Thus,  f'(t)  always 
has  one  maximum  and  one  minimum.  Clearly,  then,  f'(t)  will 
have  three  distinct  real  roots  if  and  only  if  it  is  positive 
at  t and  negative  at  t+. 


B-5 


By  substituting  from  (B-6)  into  (B-7)  it  will  be  found 


that 


f'(t  ) 


/4a  3/2 

2 + a-,  (a  maximum) 

3 vr 


and 


f' (t+) 


-4a 


3/2 


bVt 


+ (a  minimum) . 


From  these  expressions  it  can  easily  be  verified  that  when 

a1  < 0 

the  minimum  of  f'(t)  is  positive  if  and  only  if 


27ax2  > 8a23. 


(B-9 ) 


It  can  also  be  verified  that  when 


a^  > 0 

the  maximum  of  f * ( t ) is  negative  if  and  only  if  (B-9)  holds. 
Conversely,  f'(t)  will  have  a negative  minimum  and  a positive 
maximum  if  and  only  if 

8a23  > 27a:2.  (B-10) 

Since  f'(t)  can  vanish  at  three  distinct  values  of  t only 
when  its  minimum  is  negative  and  its  maximum  is  positive  (B-10) 
is  a necessary  and  sufficient  condition  that  (B-8)  have  three 
distinct  roots.  As  remarked  earlier,  this  is  also  a necessary 
condition  that  (B-4)  have  more  than  one  positive  root. 

D.  NECESSARY  CONDITION  FOR  A HAZARD 

While  (B-10)  is  a necessary  condition  that  (B-4)  have 
more  than  one  positive  root,  it  is  not  sufficient.  Moreover, 
assuming  that  (B-10)  is  satisfied,  not  only  must  additional 
conditions  be  met  before  it  can  be  concluded  that  more  than 


B-6 


I 


I 


i 


one  positive  root  of  (B-4)  exists,  but  even  when  that  fact  is 
established  the  possibility  of  a collision  before  the  time  t 
will  still  remain  open.  In  order  to  assert  that  such  a colli- 
sion is  possible,  so  that  a CAS  alarm  should  occur,  at  least 
one  of  the  positive  roots  of  (B-4)  must  be  known  to  lie  between 

0 and  t . 
e 

Before  proceeding  with  a more  detailed  analysis  of  the 
biquadratic  polynomial  f(t)  it  will  be  convenient  to  change 
variables  in  ord^r  to  reduce  the  number  of  its  coefficients 
from  three  to  two.  This  purpose  is  achieved  by  substituting 

t = s 

into  f(t)  and  dividing  the  result  by  a^ . It  may  be  observed 
that  when  (B-10)  is  satisfied  as  assumed  cannot  be  zero; 
therefore,  this  division  is  possible.  The  resulting  polynomial 
g(s)  will  have  the  form 

g ( s ) = s1*  - s2  + a^s  - aQ  = 0,  (B-ll) 


where 


In  terms  of  the  new  coefficients  (B-10)  can  be  expressed  In  the 
form 


2 


< 


8 

TT  ’ 


or 

A1  < 0.5^3.  ( B-12  ) 

Since  g(s)  Is  negative  at  s=o  and  becomes  arbitrarily  large 
for  large  s,  if  the  equation 

g(s)  = o 

has  more  than  one  positive  root,  g(s)  must  have  a maximum  at 


w 
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some  positive  value  of  s.  The  maximum  occurs  when  the 
second  derivative  g"(s)  is  negative,  i.e.,  when 

12sm2  - 2 < 0, 

or 

0 < sM  < _L=  • (B-13) 

•y6 

Also,  of  course, 

g'(sM)  = 4SM3  ’ 2SM  + A1  = 0;  1-e-* 

= 2S|vj  — (B— 14) 

From  (B-ll)  and  (B-12)  it  follows  that  a positive  maximum  can 
occur  only  if  A^  satisfies 

° < A1  < , (B-15) 

a condition  already  partly  known  from  (B-12). 

Even  when  g(s)  has  a maximum  and  is  positive,  (B-ll) 
will  have  more  than  one  positive  root  if  and  only  if  g(sM)  is 
non  negative.  After  substituting  for  A^  from  (B-14)  into  g(s^) 
this  condition  becomes 

4 2 

g(sM)  = -3sm  + sm  - Aq  i 0, 


or 


A_  s Sjyj  -3sm  = h(sM) 


3M 


5M ' 


(B-16) 


For  the  whole  possible  range  of  values  for  s^  given  by 
(B-13)  the  function  h(sM)  has  a positive  derivative  and  there- 
fore is  always  increasing.  Thus  if  (B-16)  is  satisfied  for  any 


sM  it  must  also  be  satisfied  when 

X • 6 • ) 


Ao  < F 


3 

3F 


1 

F 

1 

12 


is  substituted  for  s, 


M ’ 

(B-17) 


Let  s be  the  smallest  positive  root  of  the  equation 
o 

Ao  = h(so)  = sQ2  - 3so4; 


(B-18) 


B-8 


(B-19) 


! 


1 


1 - J 1-12A 

V— ^ 

which  must  be  real  because  of  (B-17).  Since  h(s)  is  increas- 
ing over  the  range  of  values  (B-l6)  will  be  satisfied  if 
and  only  if 

so  * SM> 


1 . e . , if  and  only  if 


2s, 


* 2s 


M 


Therefore,  g(s)  has  more  than  one  positive  root  if  and  only  if 

A1  * 2sq  - 4s03,  (B-20) 

where  sq  is  given  by  (B-19). 

When  (B-20)  is  satisfied  the  maximum  of  g(s)  is  either 
positive  or  zero.  In  the  first  case  (B-ll)  has  three  positive 
roots  and  in  the  second  case  (B-ll)  has  two  positive  roots, 
the  smaller  of  which  occurs  at  s^  where  the  maximum  also  occurs. 

If  t is  less  than  the  first  positive  root  of  (B-4)  no 
0 

alarm  should  occur.  In  this  case,  since  f(0)  is  negative  and 
the  maximum  of  f(t)  occurs  after  the  first  root  of  (B-4),  the 
slope  of  f(t)  at  t is  positive;  i.e., 

f'(t  ) = 4tg3  - 2a2t0  + a1  > 0.  (B-21) 


In  addition,  since  the  maximum  of  f(t)  occurs  before  the  posi- 
tive inflection  point  the  second  derivative  of  f(t)  does  not 
change  sign  between  zero  and  the  maximum  and  therefore  would 
have  to  be  negative  at  tg;  i.e., 

f"(t  ) = 12te2  - 2a2  < 0.  ( B-22 ) 

When  (B-4)  has  more  than  one  positive  root,  i.e.,  (B-ll)  has 
more  than  one  positive  root  and  either  of  the  conditions  (B-21) 


B-9 


or  (B-22)  are  violated,  (B-4)  will  have  a root  between  zero 
and  t and  therefore  an  alarm  should  occur. 
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